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We discuss the interdependence of resource state, measurement setting and temporal order in mea- 
surement-based quantum computation. The possible temporal orders of measurement events are 
constrained by the principle that the randomness inherent in quantum measurement should not affect 
the outcome of the computation. We provide a classification for all temporal relations among mea- 
surement events compatible with a given initial stabilizer state and measurement setting, in terms of a 
matroid. Conversely, we show that classical processing relations necessary for turning the local mea- 
surement outcomes into computational output determine the resource state and measurement setting 
up to local equivalence. Further, we find a symmetry transformation related to local complementation 
that leaves the temporal relations invariant. 

1 Introduction 

Quantum states can change over time in two fundamentally different ways, unitary evolution and mea- 
surement. The former is deterministic and reversible whereas the latter is probabilistic and irreversible. 
Quantum computation 1 1 1 can be built on either. Employing unitary evolution leads to the quantum 
circuit model ||2]|3l (the standard model of quantum computation), and is also present in adiabatic quan- 
tum computation [4] . Using measurement as the central tool to drive the computation leads to various 
measurement and teleportation-based schemes fHIHlTJ. 

Here we consider the one-way quantum computer (MBQC) [6|, a measurement based scheme of 
universal quantum computation. Therein, the process of computation is driven by local measurements 
on an entangled state of many qubits. Since the entanglement monotonically decreases over the course of 
computation, the initial state can be viewed as an entanglement resource. Known examples of universal 
resources include certain graph states |8| and ground states of two-body Hamiltonians [9], among them 
two-dimensional AKLT states of spin 3/2 [10, 11]. A classification of universal resource states is to 
date not available, but it is known that only a tiny fraction of quantum states can possibly be universal 

iniiii. 

The object of study in this paper is the temporal ordering of measurements in MBQC, and what de- 
termines it. We point out explicitly that here we do not consider a scenario where a specific quantum 
algorithm is given in the circuit model and translated into MBQC. In that case, the temporal order of 
measurements can be straightforwardly obtained from the quantum circuit by a known mapping. Rather, 
we look at the initial quantum state and the planes of the Bloch sphere in which the locally measured ob- 
servables reside, with no further information provided. We ask which temporal orders of measurements 
are compatible with this information. It may at first sight appear surprising that there is any constraint at 
all. However, the resource state and temporal order in MBQC become interrelated through the principle 
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that randomness of measurement outcomes must be prevented from affecting the logical processing. If the 
resource is a stabilizer state then all constraints on temporal order follow from the symmetries inferred 
by its stabilizer group. The purpose of this paper is to work out the consequences of this connection. 

In addition to the theory of MBQC itself, we may investigate MBQC temporal order towards a 
different goal: toy models for generating temporal order from none. Consider the case where the resource 
is a graph state | G) . How much of the temporal order follows from the interaction graph G? The graph G 
is an undirected object whereas the partial order is directed. Thus, if G constrained the temporal orders 
severely, MBQC would provide a mechanism for generating temporal order. A quantum theory of gravity 
must achieve this in a far more complicated setting. 

Previous work on MBQC temporal order has shown that for resource graph states \G) the temporal 
order of measurements is fully specified if the sets / and O of first- and last-measurable qubits are known 
lfl4l . ifTS^ l. But not every pair /, O leads to a possible temporal order. A condition on admissible pairs /, 
O in terms of the adjacency matrix of G has been given 116]. 

In this paper, we present the following results, (i) We show that all transitive temporal relations which 
prevent measurement randomness from affecting the logical processing correspond to bases of a matroid 
derived from a resource stabilizer state and the local measurement bases, (ii) It is known that the adaption 
of measurement bases according to previously obtained measurement outcomes and the extraction of the 
computational result are governed by linear- processing relations. Further, in all known schemes, the 
Bloch vectors specifying the measured local observables lie in the {X,Y), iY,Z) or (Z,X) equators of the 
Bloch sphere which are called "measurement planes". Here we show that the linear processing relations 
determine the resource stabilizer state and the measurement planes up to equivalence. For this result to 
hold we need to slightly extend the previously known processing relations by introducing gauge variable 
that affect the processing but not the probability distribution of the computational output (See Section|3]l. 
(iii) We identify a transformation that leaves the temporal order of an MBQC unchanged and which is 
a slight generalization of local complementation. Local complementation has previously been found 
useful for the discussion of local equivalence among graph states |[T7l -|[T9l. 

The remainder of this paper is organized as follows. In Section [2] we provide the necessary back- 
ground and notation, and briefly review prior work on temporal order in MBQC. In Section |3] we in- 
troduce the aforementioned gauge degrees of freedom and a resulting normal form of the resource state 
stabilizer. Based on this normal form, in Section[4]we derive our two main results on the interdependence 
of resource state and temporal order, (i) and (ii). In Section |5] we introduce a symmetry transformation 
that leaves the temporal relation in any given MBQC invariant and which is related to local complemen- 
tation. In Section [6] we show that MBQC mimics certain aspects of General Relativity f20l l?HI. We 
present an MBQC analogue of Malament's theorem |22| and discuss the emergence of an event horizon. 
In Section [7] we conclude and point out open questions. 



2 Background 

In this section we review some basic facts about measurement-based quantum computation, essential 
definitions for the discussion of MBQC temporal order, as well as previous work in this area. The 
scheme of MBQC itself, with a proof of its computational universality, is not reviewed here since various 
articles on this subject exist in the literature ||6]|, ||23l - ETl Also, we require familiarity with the stabilizer 
formalism 1281. 
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2.1 MBQC and cluster states 

In MBQC, the process of computation is driven by local (=l-qubit) measurements on an initial highly 
entangled state, generally taken to be a so-called cluster state. The local measurements can only reduce 
entanglement, and therefore all entanglement needed for the computation must come from the initial 
state. For this reason, the initial state is often called the 'resource stateQ 

The computational power of a given MBQC strongly depends on the choice of the initial resource 
state. For example, a local resource state has obviously no computational power. Other states may be 
used for a restricted class of computations. Two-or-higher dimensional cluster states of unbounded size 
have the property that they enable universal quantum computation. That is, any quantum computation 
can be realized on such a state, by suitable choice of the local measurement bases. 

We now define graph states, to be used later on, and cluster states as a subclass thereof. 

Definition 1 (Graph states and cluster states.). Be G a graph with vertex set V{G) and edge set E{G), 
such that there is one qubit for each vertex a G V{G). Then, the graph state \G) is the unique (up to 
global phase) joint eigenstate, \G) = Ka \G), for all a G V{G), of the operators 

Ka = a^^ (g) af'. (1) 

b\{a,b)eE{G) 

The cluster state \^_^) is a graph state where the corresponding graph is a d-dimensional lattice 
In the standard scheme [6], the measured observables are of the form 

0«M =cos(p«ai"' + (-l)^''sin(p,aj''\ (2) 

with a the qubit to which the measurement is applied, and qa G depending on outcomes of (earlier) 
measurements on other qubit locations. 

2.2 Temporal order in MBQC 

As noted above, the temporal order of measurement in MBQC is a consequence of the randomness of 
measurement outcomes. By adjusting measurement bases according to measurement outcomes obtained 
on other qubits, this randomness inherent in quantum mechanical measurement can be kept from creep- 
ing into the logical processing Q. If the measurement outcome of qubit a influences the choice of 
measurement basis for qubit h, clearly, qubit a must be measured before qubit h can. This is how a 
temporal order among the measurement events arises. 

We now generalize the above scenario of measuring observables of form Eq. Q on cluster states. 
Namely, we now consider general stabilizer states as resources, with supp(|*F)) = Q. and stabilizer 
group oS^(|*F)). Furthermore, we generalize the local observables whose measurement drives the com- 
putation from Eq. (|2]) to 

Oa [<?J = cos 9, C7^"^ + (- 1 sin 9« o^^ , Va e (3) 

with a0,Gs(j, (J(j> £ {X,Y,Z}. Therein, the measurement angles % are in the range —n/2 < % < 7l/2, 
and qa G 1^2 may depend on measurement outcomes from several (other) qubits in Q.. 

'Large entanglement is a necessary I29II30I but not sufficient condition for the usefulness of a quantum state as resource in 
MBQC. Somewhat paradoxically, quantum states exist that are too entangled to be useful 1121 . 
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Definition 2 (Measurement plane). For every qubit a G supp{\^)), the measurement plane at a is the 
ordered pair [o^'\ol^^]. 

We define a third Pauli operator, Os = iOsip <70 ■ As we will see shortly, the Pauli operators and Os 
are useful because of the relations 

a^O[q]al = 0[qe\], a,0[<i]<^I = -0[q]. (4) 

The basic mechanism of accounting for an "undesired" measurement outcome is the following. Sup- 
pose on some qubit a G Q., instead of the "desired" post-measurement state \ (Pa)a the "undesired" post- 
measurement state \(Pa)a has been obtained. The goal is to get the computation back on track by only 
adjusting the subsequent measurements. To do that, we require a stabilizer operator K{a) S 
with the following properties [14]: (1) K{a) has support only on a and the yet unmeasured qubits, and 

(2) K{a)\a = cii"'. Recall that (Js\(pa) = \(Pa) for the eigenstates \(pa{qa,Sa)) of the local measured ob- 
servable 0[^a], c.f. Eq. (Q. Denote by J^{a) and J^{a) the past and future of a, respectively. Then, 



(5) 



Therein, the first equality follows from K{a) £ and the second from the above properties (1) 

and (2). 

Since the overlaps between local states (representing the local measurements) with the resource 
state contain all information about the computation, we thus find that we can correct for "unde- 
sired" outcomes by (a) adjusting measurement bases of future measurements (caused by tensor factors 
a,j, in K{a) \ .^(a)) and (b) re-interpretation of measurement outcomes (caused by tensor factors Og in 

Example: Consider MBQC on a cluster state IO3) of three qubits on a line, each measured in the 
[cJ^, ay]-plane. That is, for all three qubits = Z and = X. (Here and from now on, we use the 
shorthand X = a^, Y = Oy and Z = a^.) The stabilizer generators of IO3) are 

^1 = Xi0Z20h = cy^'^®cyP®/(3), 

K2 = = aP^al^'^0a^^\ (6) 

K3 = h0Z20X3 = /(i)®af^0c7^^\ 

When the three cluster qubits are measured in the order 1 -< 2 ^ 3, the corresponding quantum circuit is 

m 



\+>-a — u,—u^\x (7) 



Using Eq. (|5]l, here we show that if qubits 1 , 2, 3 are measured in the order 1 -< 2 -< 3, then the randomness 
of the measurement outcomes on qubits 1 and 2 can be corrected for. First we consider the stabilizer 
operator K2 = Os^'^ (8> o^'' ai^'' =: K{1). If inserted in the state overlap of Eq. (5 1, the measurement 
outcome of qubits 1 and 3 is flipped, as well as the measurement basis at qubit 2. Since qubits 2 and 3 
are yet unmeasured when qubit 1 is measured, this is a valid correction operation for qubit 1 ; hence the 
notation K{1). Similarly, = 7^^) ai^^ G^'' =: K{2) can be used as correction operation for qubit 
2. K{2) flips the measurement outcome of qubit 2 and the measurement basis of qubit 3. Since qubit 3 is 
yet unmeasured when qubit 2 is measured, this corresponds to a valid correction operation. 
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The above argument also works in reverse. If the correction operations ^(1) and ^(2) are used, then 
1 -< 2 -< 3 follows. K{1) implies that the measurement basis of qubit 2 depends on the measurement out- 
come at qubit 1, hence 1 -< 2. Note that 1^3 does not yet follow! K{1) does not affect the measurement 
basis at qubit 3. Only the meaning of the eigenstates is interchanged, which by itself does not require 
qubit 3 to be measured after qubit 1 . The interpretation of the measurement outcome may take place long 
after the measurement itself has taken place. 

K3 = K{2) implies that the measurement basis of qubit 3 depends on the measurement outcome of 
qubit 2, and hence 2 -< 3. Both relations taken together yield 1 -< 2 ^ 3. 

From the equivalence with the circuit of Eq. (j7]) one would expect one bit of classical output. Indeed, 
if no correction operations need to be used, the eigenvalue measured at the output of the circuit corre- 
sponds to the eigenvalue A3 measured on qubit 3 of the cluster. Now recall that K{\), applied conditioned 
upon Ai = — 1, flips A3. Therefore, with or without corrections, the eigenvalue measured in the circuit 
Eq. ^ equals A1A3. Or, in binary notation Ai = (—1)'' , A3 = (— 1)^\ the single bit of classical output 
takes the value i'l + ^'3 mod 2. 

Note that we have not made any use of K\ = o^^^ (8) cii^' in the above argument. Still, K\ has a role 
to play, as we discuss in Section|3] 

2.3 Influence matrix, forward and backward cones 

To counteract the randomness of measurement outcomes, two measurement settings (i.e., bases) per qubit 
suffice, which may be labeled by = and = 1, respectively, for each qubit a. The measurement 
settings may collectively be described by a binary vector q, [q]^ = qa for all a^Q., and the measurement 
outcomes by a binary vector s, [s]a = Sa, forall a G It turns out that the relation between measurement 
bases q and measurement outcomes s is linear OTl . 

q = rs mod 2, (8) 

with T a binary matrix. We call T the influence matrix. 

The set of all qubits b whose measurement basis must be adjusted according to the measurement 
outcome on a is denoted as ihe. forward cone of a qubit a. Similarly, the backward cone of a qubit b is 
the set of all those qubits a whose measurement outcome influence the measurement basis at b. More 
formally, with Eq. ([8]), 

Definition 3 (Forward and backward cones). For any a £ Q. the forward conefc{a) is given by 

fc{a):={bea\dqt,/dsa=l]. (9) 

For any b G D.the backward cone bc{b) is given by 

bc{b):={a(^Q\dqh/dsa = l}. (10) 

We denote the characteristic vectors of fc{a) and bc{b) by fc(a) and hc{b), respectively. Then, the 
influence matrix T takes the form 
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The influence matrix T generates a temporal relation among the measurement events under transitivity. 
We say a <b {a precedes ft) if G fc{a). A priori, it is not forbidden that for two qubits a,b ^ D., a ^ b 
and b ^ a. However, such a computation could not be run deterministically in a world like ours where 
time progresses linearly. An MBQC is deterministically runnable if the temporal relation between 
the measurement events is a strict partial order. 

Definition 4 (Strict partial order). A strict partial order is a relation among the elements a £ D. with the 
following properties 

a ^a, ^aeil, (irreflexivity) 

a^b^^by^a, \/a,b £ (antisymmetry) (12) 

a <b,b < c a < c, \/a,b^c£Q.. (transitivity) 

Definition 5 (Input and output sets.). For a given MBQC, the input set I QQ. is the set of qubits whose 
backward cones are empty, I = {a £Q\ bc{a) = 0}. The output set O D. is the set of qubits whose 
forward cones are empty, O = {b £ D.,fc{b) = 0}. 

That is, with respect to a given temporal relation among the measurement events, / is the maximal 
set of qubits which can be measured first, and O is the maximal set of qubits which can be measured last. 

Regarding the computational output, for the purpose of this paper we are exclusively interested in 
MBQCs for which the computational result is a classical bit string. That is, every qubit in Q. is measured. 
Then, as a consequence of the randomness of individual measurement outcomes, the classical output o of 
an MBQC is given by correlations among measurement outcomes. Again, the relation between classical 
output and measurement outcomes is linear, 

o = Zs mod 2, (13) 

for a suitable binary matrix Z. 

Example: To illustrate the above notions, we briefly return to the three-qubit cluster state example of 



Section 2.2 From the previous discussion we find that 





mod 2. (14) 



Therefore, /c(l) = {2}, /c(2) = {3}, /c(3) = and bc{l) = 0, bc{2) = {1}, bc{3) = {2}. Hence, 
/ = {1} and O = {3}. Also, 1^2 and 2^3. The latter two relations generate a third under transitivity, 
namely 1^3. Asymmetry and irreflexivity are obeyed in this example. Regarding the single bit of 



output in this computation, the matrix Z of Eq. (13lisZ = (101 



2.4 Brief review of prior work on MBQC temporal order 

MBQC temporal order as an (almost) emergent phenomenon. In II14I1 the following question is asked: 
"Given graph G and the set Z of measurement planes for all vertices, can the temporal order of measure- 
ments in MBQC with a graph state \G) be uniquely reconstructed from this information?" The graph G 
and the measurement planes £ are undirected objects. Thus, if the answer to this question was yes, then 
temporal order in MBQC were truly emergent. 

However, it turns out that the pair G,S does not specify the temporal order of measurements in 
MBQC uniquely; there are in general a number of consistent temporal orders respecting the requirement 
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that the randomness of measurement outcomes should not affect the logical processing. One may then 
ask how constraining on temporal order this requirement actually is. To this question, the following 
answer is provided by [MJ : If in addition to G and Z the set / of first-measurable and the set O of last- 
measurable qubits is known, then the complete temporal order (if existing) can be uniquely reconstructed 
from this information. Thus, MBQC temporal order is not emergent in the strict sense; a seed /, O must 
be provided in addition to G and Z, and the complete temporal order then follows. 

But not every pair /, O will lead to a consistent temporal order. The question that now arises is which 
pairs /, O do. For the case where the stabilizer resource state is a graph state and all qubits are measured 
in the [Z,F]-plane|^then the answer to this question is given in ||T6l . Denote by the adjacency matrix 
of the graph G describing the resource state |G), and by Ag|o' x/' the submatrix of Ac where the rows are 
restricted to O'^ := Q\0 and the columns are restricted to I'^. Then, the pair /, O leads to a partial order 
of measurement events in MBQC iff there exists a matrix T such that AgIo^xfT = I, and T is free of 
cycles (that is Taa = 0,Va, TahTba = 0,\/a,b, TatTbcTca = 0,\/a,b,c, etc). The resulting temporal order is 
generated by T under transitivity. 



a) b) 




Figure 1: Closed time-like curves in MBQC. a) Bennett, Schumacher and Svetlichny's post-selection 
model ll33l[34]| of CTCs (left: circuit with wires 'going backwards in time', right: implementation thereof 
using teleportation and post-selection), b) Nested forward cones in the MBQC equivalent of the telepor- 
tation circuit in (a). 

MBQC and closed time-like curves. In ll32l it is shown that MBQC encompasses the post-selection 
model of closed time-like curves (CTCs) proposed by Bennett, Schumacher [33] and Svetlichny ll34ll . 
The CTCs arise from circuits such as the one displayed in Fig. [T^, translated into MBQC. The result 
are forward cones with the property b G fc{a) Aa S fc{b), for two suitably chosen qubits a,b ^ Q.; see 
Fig.[T]5. Such nested forward cones are an obstruction to deterministic runnability of MBQC, but mimic 
closed time-like curves of General Relativity in the MBQC setting. 

3 Gauge degrees of freedom 

Here we introduce the notion of "gauge transformations" acting on a given quantum computation. These 
transformations exist for both the circuit model and MBQC. 

3.1 Gauge transformations in the circuit model 

To obtain an intuition for the gauge transformations introduced here, it is instructive to first inspect them 
in the circuit model. Specifically, we consider a quantum circuit which consists of (1) the preparation of 

^The former condition alone is not a restriction, since all stabilizer states are local Clifford equivalent to graph states 1351 . 
but both conditions jointly are. 
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a) b) 
/ 




Figure 2: Two symmetry transformations, a) Gauge transformation in the circuit model. For any logical 
qubit I, an identity I = ax ''Gx '' is inserted into the circuit next to the input. The left Ox is propagated 
backwards in time, and absorbed by the input state |+). The right a,- is propagated forward in time, 
flipping rotation angles and, potentially, measurement outcomes in the passing, b) Flipping a measure- 
ment plane in MBQC. In the measurement plane [00,0^0], the Pauli operator 00 is distinguished over 
(Jsip because the rule for adjusting a local observable O for measurement is 0[q = 1] = CJ0O[g = OjcrJ. 

If, for a qubit a, Taa = (with T the influence matrix) then the exchange o^'^ < — > G^^^ is a symmetry 
transformation for the given MBQC. 



the quantum register in the initial state ®JLi 1+)/. (2) unitary evolution composed of, say, CNOT gates 
and one-qubit rotations about the X- and Z-axes, and (3) local measurements for readout. Such a circuit 
is displayed in Fig. |2]above. Then, into every qubit line individually, we may insert an identity / = Oxf^x 
next to the input; See Fig. [2] The left Ox is propagated backwards in time until absorbed by the input 
state |+). The right CJv is propagated forward in time, flipping rotation angles and readout-measurement 
outcomes in the passing. 

This transformation is an equivalence transformation, since it is caused by the insertion of an identity 
gate into the circuit of Fig. [2] It changes the sign for certain rotation angles, i.e. when angles are counted 
positive or negative. Specifically, for the z-rotation gates next to each input qubit we can individually 
choose our convention for which rotation angles are called positive or negative, respectively. Once those 
signs are fixed on the input side, they are fixed throughout the circuit. Changing this reference affects the 
procedure of computation, but leaves the distribution of computational results unchanged. We therefore 
call it a gauge transformation. 



3.2 The gauge transformations in the measurement-based model 



Translating the above discussed gauge transformations from the circuit model into MBQC it is easily 
seen that the above relations Eq. ([8]) and ( 13 1 are incomplete. We find the more general relations 



q = rs + //g mod 2, 
o = Zs + /?g mod 2, 



(15a) 
(15b) 



with g a choice of gauge. 

Under a change of g, the measurement bases in a particular MBQC change, but the probability 
distribution for the classical output values remains unchanged. Of course, knowledge of the g-dependent 
extra parts in Eq. ( 15a|15b i is not necessary to run any given MBQC, since g = is always a valid choice. 
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However, the presence of the extra terms in the processing relations strengthens their interdependence 
with the resource state, which is the reason why we discuss them here. 

We may now want to derive the generalized processing relations Eq. ( 15a|15b ) directly in MBQC, 



without reference to the circuit model. Before discussing the general case, we return to the specific 



example of the three-qubit cluster state in Section 2.2 



Example: Consider the product K\K3 = al'"* (g) aP^ of ^i, ^3 in Eq. (|6|. When used in Eq. (|5|, the 



effect of this stabilizer element is to flip the measurement bases of qubits 



[ and 3. Hence, the relation 



Eq. ( [T4| ) generalizes to 

f qi \ f \ f si \ f \ \ 

mod 2. (16) 





This is precisely what one would expect from the insertion of GxGx next to the input of the equivalent 
quantum circuit in Eq. (|7]l. 

We now turn to the general case. Via Eq. ([Sjl, the stabilizer group acts on s and q. Denote the 

post-measurement state of qubit a G 11 by \sa,qa)a, with Sa the measurement outcome and qa specifying 
the measurement basis. Then, as in Eq. ([5]l, 

)a{Sa,qa\]\'i')= 1 1^ a{Sa,qa\] Kl^') = [ , ^„ |^ ) , V/:G^(|»F)). (17) 
\aeCl J \aeQ. J \aeQ. J 

Since Gs\s,q) = \s(Bl,q) and a^\s,q) = \s,q(Bl), for a stabilizer element A" = <S)aeiii^s''^y''' {(y^"^)^'" the 
action of Gk on s, q is 

q — ;> q + w mod 2. ^ ^ 

Again, nothing changes by the insertion of a stabilizer (identity) operator into the state overlap of Eq. Q, 
and transformations Gk are therefore equivalence transformations. They can be used to constrain the 
possible temporal orders in MBQC, as we discuss explicitly in Appendix [A] 



3.3 Closed time-like curves - good or bad? 

Partial temporal orders for an MBQC with a given resource stabilizer state and fixed measurement planes 
are the solution to two constraints, namely 

1 . Every qubit in the resource state must have a forward con^ 

2. The forward cones generate an irreflexive and antisymmetric relation under transitivity. 

In this paper we allow closed time like curves (CTC) in the computation; we are interested in classifying 
all transitive temporal relations consistent with a given resource stabilizer state and set of measurement 
planes. We do not restrict to partial orders per se, and therefore drop the above Condition [2] The 
remaining Condition [T] may, at first sight, hardly seem to pose any constraints at all. However, it does. 
It imposes a self-consistency condition on each forward cone. As we make explicit in Section [6]2j this 
self-consistency condition in certain cases takes the form of a wave equation. 

A temporal relation among measurements which is not a partial order contains closed time-like 
curves, and can be implemented in familiar linear time only by employing postselection. On the other 

■^Forward cones are allowed to be empty. 
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hand, temporal relations which contain closed time-like curves have recently been found of indepen- 
dent interest [32|. They are the translation of quantum circuits with post-selection CTCs |[33l . [34] into 
MBQC. We may compare with the theory of General Relativity, where certain solutions of the Einstein 
equations contain closed time-like curves. Such solutions give rise to a host of paradoxes, and whether 
they are physical is under debate |36]. But Einstein's field equations are not abandoned because they 
allow for CTCs. 



3.4 Correction and gauge operations in the stabilizer formalism 

To state and prove our results on MBQC temporal order, we need to make a few more definitions. It turns 



out that the possible temporal relations, and indeed the classical processing relations Eq. (15ai, ( 15b i 



can be parametrized by two subsets of D., namely the computational output set Ocomp and the gauge input 
set /gauge- We define these sets next. 

We say that the measurement outcome Sa of qubit a is corrected in a given MBQC, if by insertion 



of a suitable stabilizer operator K in Eq. ( 17 1 equivalence of 5q = 1 with the reference outcome Sa = 
is established, at the cost of adjustment of measurement bases and/or re-interpretation of measurement 
outcomes on other qubits. 

We observe that if for a given MBQC all measurement outcomes Sa, a e Q., can be corrected then 
these measurement outcomes contain no information and no linear combination of them is worth out- 
putting. Thus, in general there will be a set of qubits whose measurement outcomes are not corrected. 

Definition 6 (Computational output set). For a given MBQC, the computational output set Ocomp Q ^ is 
the set of qubits whose measurement outcomes are not corrected. 

The correction operations for the qubits in a G (Ocomp)'^> c.f. Eq. (j5|l, are each implemented by 
correction operators K{a) G For each a G (C?comp)'^> ^(^) has the property that it flips the 

measurement outcome Sa at qubit a, but does not flip the measurement outcome on any other qubit in 
(Ocomp)'^- In this way, it is ensured that the correction operation is for qubit a individually. 
Definition 7 (Correction operator). For an MBQC with a given stabilizer resource state fixed set 
Z of measurement planes and computational output set Ocomp, for each a G {OcompY the corresponding 
correction operator K (a) ^ is a Pauli operator satisfying the conditions 

' ' (19) 



K{a)U G {ai^Ui;'}, 



The correction operator K{a) can be used to correct an "undesired" measurement outcome at qubit a, 
c.f. Eq. (5 1. If, for a given operator K{a), K{a)\b G {a^*\ a|^^} then, by Eq. (5 1, the measurement basis 
at qubit odepends on the measurement outcome for qubit a. In terms of the influence matrix T , 

K{a)\b G {af , ajj'} ^ [r],,, = 1, Va G (Ocomp)' , ^ G (20) 



Note that in the first line of Eq. (19 1 we allow K{a)\a = cjj^' only because we are admitting closed 

time-like curves in the present discussion. K{a)\a = cjj^^ means that the measurement basis for qubit a 
depends on the outcome Sa of the measurement of qubit a. This amounts to a closed time-like curve only 
involving qubit a (self-loop) and is an obstacle to deterministic runnability. 

Because the qubits a G Ocomp have no correction operations, Tba = for all b ^ Q., and thus 

Ocomp Q O. (21) 



R. Raussendorf, R Sarvepalli, T.-C. Wei, and P. Haghnegahdar 



229 



Thus, Ocomp and the correction operators {K{a) , a G (Ocomp)'^} completely determine the influence matrix 
T and hence the temporal relation among the measurements. 



Conversely, if Ocomp and {K{a),a € (Ocomp)'^} are unknown, then the constraints Eq. ( 19 1 pose self- 
consistency conditions on them. We discuss these conditions further below. 

As we have seen in the concrete three-qubit example above, for a given MBQC the relation between 
the choice q of measurement bases and the measurement outcomes s in general allows for an offset term 



Hg, c.f. Eq. ( 15a I. Thus, the measurement bases for a certain set of qubits can be freely chosen until the 
initially arbitrary g becomes fixed. This observation leads to 

Definition 8 (Gauge input set). For a given MBQC with input set I, the gauge input set Igauge ^1 is a set 
of qubits such that for each i G Igauge the parameter qi specifying the locally measured observable Oi[qi\ 
can be freely chosen. 

Analogously to the correction operations, there will be gauge operators which implement the 'cor- 
rections' of measurement bases for the qubits in /gauge- The definition of the gauge operators ensures that 
for all / G /gauge the corresponding g,- can be changed individually without changing the others. 

Definition 9 (Gauge operators). For an MBQC with a given stabilizer resource state fixed set £ of 
measurement planes, computational output set Ocomp and gauge input set Igauge, for each i G Igauge, the 
corresponding gauge operator K(i) G o5^(|*F)) is a Pauli operator satisfying the conditions 

m\i = a«, 

^(01; = yj€{Igaugen{0,o,npr)\i, 
K{i)\, G {/«,af}, yk€{{IgaugeYn{0,o,npr)\i, 

K{i)\, G {/('), ai')}, yie{IgaugenOcomp)\i. 



Like previously for the correction operations, Eq. ( [22] ) poses self-consistency condition on the possi- 
ble sets /gauge and {K{i),i G /gauge}- 



Lemma 1. The correction operators K{a) of Eq. (191, a G {OcompY, and the gauge operators K{i) of 
Eq. {22 ), i G Igauge, are independent. That is, for all sets J C Igauge, L C {OcompY with 7 / V L 7^ 0, 



KiJ,L) := UiejKii) UaeiKia) / /(^). 

Proof of Lemma [7] (indirect) Assume that there exists a pair J, L, with J ^'^M L such that 
k{J,L) = /("). Then, for all b G (Ocomp)S K{J,L)\t = 1. Now, with Eq. ([l9i, K{b)\t G {ai^', aj^^^}. 



and K{c)\b G {aj^',/(*)} for all c G (Ocomp)'\^- With Eq. \22\, K{i)\b G {a)''\/W} for all / G I, 



'gauge - 



Therefore, no other K{-), K{-) can cancel a a, -contribution from K{b) to K{J,L). Hence, b ^Lior all 
b G (Ocomp)'^, and thus L = d). By an analogous argument, no K{-) can cancel the (7^'^ -contribution of 
K{i) to K{J,L), hence / J for all / G /gauge, and 7 = 0. Thus, K{J,L) = /("^ =^ J,L = (d. Contradiction. 
Hence, the K{a), K{i) are independent. □ 

Remark: For any given MBQC with fixed temporal relation, the computational output set Ocomp is a 



subset of the output set O, see Eq. (21 1. But Ocomp and O are not necessarily equal. Likewise, /gauge ^ / 
by definition, but /gauge and / may not be equal. To illustrate this point, we consider the following two 
examples. 



Example 1. Consider the 3-qubit cluster state of Section 2.2 with measurement planes [X/Y] for 
all three qubits. We consider the correction operators ^(1) = K2, K{2) = Kt, for qubits 1 and 2, and 
gauge operator ^(1) = KiKi,. Ocomp = {3} and /gauge = {1} are then permitted by Eqs. (19 1 and (22 1, 



230 



Symmetry constraints on temporal order in measurement-based quantum computation 



respectively. As was discussed previously for the above choice of correction operations, / = { 1 } and 
O = {3}. Thus, in the present example /gauge = / and Ocomp = O. 

Example 2. Consider a Greenberger-Horne-Zeilinger state \GHZ) = (|000) + 1 11 1))/ \/2 as resource 
state, with all three qubits measured in a basis in the [X,F] -plane. We use the stabilizer elements K{\) := 
Z1Z3 = ai^'aP' and ^(2) := Z2Z3 = ai^^aP^ as correction operations for qubits 1 and 2, and A'(l) := 

X1X2X3 = o^^'' o^^^ as gauge operator. Then, the choice /gauge = {1} and Ocomp = {3} is admitted by 
Eqs. ( 19 1 and ( [22] ), respectively. On the other hand, this is an example of a temporarily flat MQC, T = 0. 



Therefore, / = 6> = {1,2,3}, and Ocomp / O, /gauge / /■ 
Lemma 2. For any MBQC on a stabilizer state, \Igauge\ ^ \Ocomp\ 
We prove Lemma[2]in Section 3.5 



Definition 10. A pair 1 gauge, Ocmp is called extremal iff\Igauge\ = \Ocoinp\- 

As will become clear in the next section, extremal pairs /gauge > Ocomp are easier to handle than general 
pairs, and are not very restrictive (c.f. Theorem[2]). 



We still need to relate the classical output vector o appearing in Eq. ( 15b 1 to the set Ocomp- To this 
end, we make the following 

Definition 11 (Optimal classical output). A classical output vector o, with processing relations o = 
Zs+Rg, is optimal iff the following conditions hold 

1. Maximality: Upon left-multiplication by an invertible matrix, Z can be brought into a unique 
normal form 

Z~(Z|/), (23) 
where the column split is between {OcompYlOcomp, cind 

2. Determinism: For the matrix Z in Eq. pi]), 



[Zhy=l^/^(j)|,-GK\cTW}. (24) 
In other words Eq. ([24| informs us that the jth column of Z is simply the restriction of the support of 

K{i) to Ocomp- 

The reason for defining an 'optimal classical output' besides a 'classical output' is the following: 
One could, in principle, run an MBQC perfectly deterministically and then choose o such that nothing 
is outputted at all, or all outputted bits, independent of the measurement angles chosen, are zero guar- 
anteed or perfectly random guaranteed. The above definition of an optimal classical output eliminates 
such choices. Maximality says that there is one bit of optimal output per qubit of Ocomp- The deter- 
minism condition can be understood from the correction procedure explained in Eq. ([5]l. For a qubit 
i G (Ocomp)'^, we account for the undesired outcome Sj = 1 by inserting K{i) into the state overlap 



(<?'ioc|*I') = {fP{oc\K{i)\^) . Now consider a qubit / G Ocomp- By Eq. (23 1, Si contributes to the output 
bit Oi, the /-th bit of o. If K{i)\i G {Os^\o\'^}, then the insertion of K{j) into the overlap flips Si, i.e.. 
Si — )■ Si@\. This needs to be taken into account when reading out Si. The linear combination Sj@Si 
remains unaffected by the correction for Sj. 

We have now made the definitions needed to state and prove our results on temporal order in MBQC. 



In Section 3.5 we establish a normal form for the stabilizer generator matrix of the resource state 
This normal form is the basis for our results on MBQC temporal order, which are stated in Sections [4] 
andlAl 
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3.5 Normal form of the resource state stabilizer 



Let us briefly review which pieces of information specify an MBQC. A priori, there are four: the set of 
measurement angles, the set £ of measurement planes, the resource state and the classical process- 
ing relations Eq. ( 15a|15b i. The measurement angles entirely drop out of all our considerations about 
temporal order. Next, we observe that when specifying the measurement planes and the resource state 
separately, we really specify too much. Starting from a given pair l*!'),!^ of resource state and set of 
measurement planes, for any local Clifford unitary U, the pair U\^),U{L) obtained by applying U to 
both the measurement planes £ and the stabilizer state is again a valid pair, i.e., it consists of a set 
of a stabilizer state and a set of measurement planes. Furthermore, it amounts to exactly the same com- 
putation as the original pair. The pair is thus redundant. To remove this redundancy, we combine 
the measurement planes and the stabilizer state into the stabilizer generator matrix in the 
Otj, /ai-stabilizer basis, 

^{m = ms). (25) 



Therein, the columns to the left (right) form the G,j, (dr) part of the stabilizer generator matrix. §f 
comprises all information from the resource state and the set of measurement planes £ relevant for 
the discussion of MBQC. We do not need to know the state and the measurement planes separately. 

Now note that the correction operators K{a), a G (Ocomp)'. and the gauge operators K{i), i G /gauge, 
are all elements of the stabilizer ^{\^')), and, by Lemma[T| are independent. Thus, they either form or 
can be completed to a set of generators for o5^(|*F)). This observation leads us to the following 

Lemma 3. For any MBQC on a stabilizer state with extremal 1 gauge, Ocomp< the generator matrix $f 
o/^(|*F)) can be written in the normal form 

<y<j) CTs 

^ gauge gauge) {Ocomp) Ocomp 



( 




/ 




V ' 










(26) 



The matrices H, R, T, Z are related to the matrices H, R, T, Z governing the classical processing of 
measurement outcomes in MBQC, q = Ts + Hg mod 2 and o = Zs+ Rg mod 2, via 








T 






H 



Z= (Z|/) ,R = R. 



(27) 



Proof. By Definition [Tj a correction operator K{a) exists for every a G (Ocomp)' - By Eq. (19 1 these 
operators have no support in / I) /gauge and are 1(0, 
( I A II / I B ), where the part is split as /gauge, /i 
The gauge operators are, by definition, of the form: ( 

the correction operators. Since there are \Q.\ generators for the stabilizer of which | (Ocomp 

already accounted for, there can be at most |Ocomp| independent gauge operators, i.e., |/gauge| < |Ocomp| 
which proves Lemma |2] In the present setting, the pair /gauge, Ocomp is extremal by assumption, thus 



comp; I in number and must take the following form: 
gauge while the a,-part is spht as Ocomp, (Ocomp)''- 
/ I C II I D ) , where the column splits are as for 

are 



232 



Symmetry constraints on temporal order in measurement-based quantum computation 



these two sets of generators exhaust the stabilizer generators and we can write the stabilizer as 

-^gauge (^gauge) (f^comp) f^comp 







A 


I 


B 




[ ' 


C 





D 



(28) 



gauge 



c/. 



for suitable matrices A, B, C and D. We now need to identify these matrices. By definition, /j 
Measurement outcomes on qubits a € Ocomp are not corrected, hence fc{a) = for all a G C?comp> and 
Ocomp ^ O follows from the definition of the output set O. Then, the influence matrix T takes the form 



(29) 








T 






J gauge y 



Now consider the correction 



where the column split is (C?comp)'^| Ocomp and the row split is /gaugel (4 
operator K{a) for a € (Ocomp)'^ in the upper part of ^^(|*F)) in Eq. (28 1. We already know from Eq. ( |20| ), 



that the part of K{a) is the forward cone of a. Therefore we must have A = T . Further comparing, 
with Eq. ( [24] ) which states that the restriction of the a, part of the correction operator K{a) to Ocomp is 
the ath column of Z. But this is precisely the ath row of B, thus we infer that B = tJ . 

Next, consider row / of the lower part of '^i}^)) in Eq. (|28]). Row / is (0, ..,0, 1,0, ..,0|c^||0|d^). 
The corresponding stabilizer operator K(i), when inserted into the overlap as in Eq. (jsj), flips 

the measurement basis at qubit / and of qubits / G (/gauge with [c]/ = 1. It further flips the measurement 
outcomes at qubits m G Ocomp with [d]m = 1. Therefore, 



R = R, 



with H = C^, and 
with R = D^. 



We thus arrive at the normal form Eq. (26 1. 



□ 



4 Interdependence of resource state and temporal order 



Let us return to our discussion from the beginning of Section 3.5 on which pieces of information are 
needed to describe an MBQC. At this stage, apart from the set of measurement angles which do not 
enter our discussion, we remain with two pieces of data specifying a given MBQC, namely and 
the processing relations Eq. (15ai, (15bi. But it doesn't stop there. As the results of |[T4ll . |[T5l , |[T6l 
show, the resource state the measurement planes £ and the influence matrix T — being part of the 
classical processing relations Eq. (15ai, (15b I — are not independent. Specifically, given the sets / of 
first-measurable and O of last-measurable qubits in addition to and £, the temporal order (generated 
by T) can be worked out completely. 

Here we prove a statement about the interdependence of and the MBQC classical processing 

relations which goes the opposite direction. Namely we show that the classical processing relations 



Eq. (15ai, (15bi uniquely specify the stabilizer generator matrix i.e., the pair £ up to 

equivalence; See Theorem [4] below. Thus, only two pieces of data are needed to specify an MBQC 
that satisfies the determinism constraints, namely the measurement angles and the classical processing 
relations for the measurement outcomes. 
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A further question is whether the temporal relations compatible with a resource state and set of 
measurement planes £ fit into a common framework. In this regard, we show that the classical processing 
relations (containing the temporal order) for MBQC with a fixed resource state and set of measurement 
planes, for extremal pairs /gauge, Ocomp, are in one-to-one correspondence with the bases of c.f. 
Theorem 121 

4.1 Results 

We now present four theorems on the mutual dependence of the resource state and the classical process- 
ing relations. 

Theorem 1. Consider an MBQC on a stabilizer state with fixed measurement planes and an ex- 
tremal pair of gauge input set 1 gauge ^nd computational output set Ocomp- Then, the relations q = Ts +Hg 
mod 2, and o = Zs +Rg mod 2 for an optimal output o are unique. 

That is, once the resource state the measurement planes and /gauge, Ocomp are fixed, there is 
no freedom left to choose the classical processing relations. They are uniquely determined by the for- 
mer. In particular, for fixed stabilizer state and measurement planes, T = T (/gauge, Ocomp), H = 

//(/gauge, Ocomp) CtC. 

A corollary of Theorem[T]is that given the measurement planes and an extremal pair /gauge, Ocomp, the 
resource state uniquely determines the influence matrix T. One may ask how restrictive a condition 
the extremality of the pair /gauge, Ocomp is. In this regard, note 

Theorem 2. Consider an MBQC on a fixed resource stabilizer state for fixed measurement planes, with 
an infiuence matrix T and input and output sets I{T), 0{T), such that no qubit a £F can be individually 
gauged with respect to I(T), 0(T). Then, there exists an extremal pair Igauge Q 1, Ocomp ^ O such that 

T — T {IgQugg^iO comp\ 

The input set I{T) and the output set 0{T) which appear in Theorem[2]are uniquely specified by T 
through Definition |5] T (/gauge, Ocomp) is uniquely specified by the pair /gauge, Ocomp through Theorem [T] 

Theorem[2]states that all temporal relations for an MBQC, subject to the extra condition on the qubits 
which can be individually gauged, arise from extremal pairs /gauge, Ocomp- By establishing Theorem [2] 
we trade the condition of the pairs /gauge, Ocomp being extremal for the condition that no qubit in F can 
be individually gauged. The latter is a more meaningful condition. Suppose a qubit a in /' could be 
individually gauged wrt /, O. Then K{a) exists. For any b with K{b)\a = CT^"\ K{b) := K{b)K{a) is a 
valid correction operator for qubit b, and K{b)\a = l'^'^^ ■ Hence, a could be removed from all forward 
cones and thereby be made a qubit in /. By imposing the extra condition in Theorem |2} we exclude 
temporal relations where certain qubits could be in the input set / but aren't. 

Theorem [T] is mute on the question of which extremal pairs /gauge, Ocomp are admissible. Theorem [3] 
below describes how much freedom remains for the choice of the classical processing relations, given 

Theorem 3. For MBQC with a fixed resource stabilizer state and fixed measurement planes, the 
classical processing relations for extremal Igauge, Ocomp, specified by the matrices H, R, T, Z and the 
sets Igauge: Ocomp, dfc in one-to-one correspondence with the bases of the matroid ?^(|*P)). 

After we have justified our restriction to extremal pairs of gauge input and computational output sets 
in Theorem[2]and have characterized the set of temporal relations compatible with a given resource state 
and set of measurement planes in Theorem |3] we now return to Theorem [TJ and show that a converse 
also holds. 
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Theorem 4. Consider an MBQC on a stabilizer state with classical processing relations q = Ts + 
Hg mod 2, o = Zs + Rg mod 2 for an optimal classical output o, such that rkH = rkZ. Then the 
classical processing relations uniquely specify the stabilizer generator matrix in the 0,^/0s- 

basis, i.e. the resource stabilizer state and set £ of measurement planes up to equivalence. 

Remark: For Theorem |4] it does not matter whether or not the classical processing relations codify a 
temporal relation which is a partial order. 

Remark: There is a constructive procedure for obtaining ?f (in the ai^/a^ -basis) from the linear 
processing relations. If the processing relation did not stem from an actual computation but rather was 
"made up", the resulting may not be a valid stabilizer generator matrix. I.e., the rows of 

may correspond to Pauh operators which do not pairwise commute. 

4.2 Proofs of Theorems IM] 

Theorem[T]is an immediate consequence of Lemma|3] 

Proof of Theorem^ Assume that / is valid input set and O is a valid output set for a given MBQC. Then, 
the stabilizer generator matrix of the resource state can be written in the a^/av-basis as 



\ 



I 


r 


0' 










I 


A 


B 


c 





D 



) 



(30) 



for some matrices A, B, C and D. The influence matrix T can be obtained as 















(31) 



with the column split between / and F, and the row split between and O. The matrices B and {B\C) do 
not necessarily have maximal row-rank. By row transformations of (B|C| |0|D) we extract the dependent 
rows, and obtain 



\ 



I 


r 


0' 










I 


A 











D' 





C" 





D" 


B'" 


C" 





D'" 



\ 



(32) 



However, note that each row in the third set of rows in the above matrix can either be interpreted as a 
correction operator K{a), with non-empty forward cone fc{a), for some a G O, or a gauge operator K{i) 
for some / G F . The former is ruled out because every a £ O must have an empty forward cone. The 
latter is ruled out by the assumption that no qubit outside the set / can be individually gauged. Therefore 
that set of rows must identically vanish. 
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Since has full row rank, so does the matrix D' appearing in Eq. (32i. We may then choose 

a set AO C O such that the columns of D' indexed by AO form a maximal independent set. We set 

Oc 



'comp '■= 0\A0. Then, by further row transformations which do not affect f , the matrix in Eq. (32i can 
be converted to 



V 



/ 


F 


0' 


AO 


f^comp 









I 
















/ 


A" 


(33) 


B'" 


c'" 








D'" ^ 





B'" has full row rank by construction. We can therefore find a set /gauge ^ / such that the columns of B'" 
indexed by /gauge form a maximal independent set. For any such set /gauge we can convert the matrix in 
Eq. ( [33] ) fully into the normal form Eq. (|26j) without affecting T . 

For any of the above choices for /j 
can be extracted as 



gauge ^ / and Ocomp ^ O, the resulting influence matrix T (/gauge, Ocomp) 









T (/gauge 5 Ocomp ) — 
















with the column split between / and F, and the row split between O'^, AO and Ocomp 
comparison of Eqs. (31 1 and (34i we verify T = T{Ig 



'gauge ) f^comp) 



(34) 

0\A0. By 

□ 



Remark: Comparing Eq.([33|) with the normal form in Eq. ( |26| ), we can write the stabilizer matrix in 
a shghtly varied form that can be useful later. 







/gauge 


A/ 


/' 


0' 


AO 


^comp 




( 






















I 





















/ 


tJ 




/ 




4 











(35) 



where Z = (Zi | Z2) and where 
Proof of Theorem^ Denote by =5 



= (h[|hD. 

the set of bases of ' 



^ and by the set of extremal classical pro- 
cessing relations of form Eq. (15ai, (15bi, specified by the triple (/gauge, Ocomp,{H,R,T,Z}). Then, the 
mapping h: ^ - 



— > S/" exists and is a bijection. (1) Existence of h: By the normal form Eq. (26) of 
^ = (<I>| |S), for a given basis B{^) the sets /gauge and Ocomp are extracted as follows. A qubit / is in /gauge 
if and only if the corresponding column of <I> appears in B{W). A qubit a is in Ocomp if and only if the 
corresponding column of S does not appear in B{W). Knowing /gauge and Ocomp, {H,R,T,Z} is uniquely 
determined via Theorem[T] (2) Surjectivity of h: By definition of "extremal". (3) Injectivity of h: given 



'gauge 

and Ocomp, B{G) 



(*J>lwl%omp) is unique. 



□ 



Proof of Theorem^ We divide the proof into the following steps: i) First, with rkH = rkZ, the process- 
ing relations stem from an extremal pair /gauge, Ocomp- We show that given an extremal pair /gauge, Ocomp, 
the matrices H, R, T and Z are uniquely determined by the classical processing relations Eq. ( 15a|15b ). ii) 
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From these matrices we can derive the corresponding normal form of the resource state uniquely. Let us 
denote this by In iii) and iv) we show that the following diagram commutes, which establishes the 
equivalence of normal forms of all extremal pairs. 



(^gaugC) Ocomp) ^ (^gauge'^c 



comp 
T,H,Z,R 



M(A) 



gauge ' comp / 

t',h',z',r' 



(36) 



In Eq. ( [36l ), M(A) is a transformation on the normal form jV , dependent on A. Because of their inde- 

^comp separately in iii) and iv) 



^gauge 'irid f^comp 



pendence, we consider the transformations /gauge 
respectively. 

i) Extracting H, R, T and Z: Given any set of inputs / and outputs O by Theorem |2] we know that there 
always exists an extremal pair (/gauge, Ocomp)> where /gauge ^ /, Ocomp ^ O. By Eq. (29 1, T is uniquely 



specified by T . By the assumption of the classical output being optimal, Z can be brought into the 



unique normal form of Eq. ( 23 1 by left multiplication with an invertible matrix and (permuting the 
columns if necessary), Z can be extracted. The matrix H may not appear in Eq. ( 15a| ) in its normal 
form Eq. ( [27] ), nonetheless for an invertible matrix A, the vector q specifying the measurement 
bases is invariant under H — > //A, g — > A^'g. By definition of /gauge> every qubit in /gauge can 
be individually gauged with respect to /gauge, Ocomp- Therefore, (up to row permutations), we can 
choose A such that 

^ ' (37) 



HK 



H 



where the row split is between /gauge (upper) and 
output o is invariant under the transformation R - 



/gauge (lower). A and H are unique. The classical 
/?A', g — > (A')^'g. However, since Eqs. ( 15ai 
and (15b I refer to g in the same basis. A' is now fixed: A' = A. Then, R = /?A, with the A of Eq. (37 1. 



ii) Assembling the normal form: Using the unique matrices H, R, T, Z, by Lemma |3] we can now 
assemble the normal form Eq. (26 1 of the stabilizer generator matrix for the resource state By 



assumption of Theorem|4j the matrices H,R,T ,Z describe a vahd computation, and the normal form 



Eq. ( 26 1 derived from them must thus yield a valid description of a quantum state. In particular, all 



Pauli operators specified by the rows of 1#({H,R,T,Z}) in Eq. (26l must commute. (The rows are 
independent by design of the normal form.) We have thus constructed a description of Since H, 
R, T, Z are unique, so is 

We now proceed to construct the stabilizer from the processing relations when /gauge and 



Ocomp are not specified. From the classical processing relation Eq. (15ai we can still extract the 



input set / and the output set O, by testing which rows and columns of T identically vanish. Then, 
the possible choices for /j 



gauge and Ocomp are limited to /gauge ^ / and Ocomp ^ O by definition. 



iii) Equivalence under A : I^ 



gauge 



4auge- order to prove this we rely on the slightly variant version 



of the normal form of ^^(|*F)) as shown in Eq. (35 1 which also includes additional detail about the 
correction operators. As noted above, the different normal forms for H can be interconverted by 
right-multiplication with an invertible matrix A (change of basis for g), i.e. H, R — > HA, RA, where 
= (/|H^). Under such a transformation, the upper part of the normal form Eq. (26 1 for 



remains unchanged, and the lower part is transformed (//^||0|/?^) — > A^(//^| |0|R7j. Invertible 
row transformations onW{\^)) leave unchanged, and is thus independent on the precise 

choice of /gauge ^ /■ 
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iv) Equivalence under A : Ocomp — ^ f^comp- Proving that a different choice of Ocomp does not change the 
stabilizer state is a little more complicated. We proceed in the following manner. From Theorem |3] 
we know that /gauge U Ocomp "^he bases of a matroid. Therefore for two distinct computational 
output sets Ocomp and Ocomp' there exists another computational output set O'^^^p = Ocomp \ {^l U { j}, 
where / G Ocomp \ 0"onip and j G Ocomp \ ^comp- Therefore, it suffices if we show that the stabilizer 
state does not change if we change the computational output set from Ocomp to Ocomp- Assume that 
the classical relation for the computational output set Ocomp is given as 

o = Zs + /?g = (Zi|Z2|/)s + Rg. (38) 

where the column split of Z is between C?comp' and Ocomp- The corresponding normal form 
(with the column split in 00 -part /gauge l^gauge ) is 








/ 
















/ 

















(39) 



where 



Suppose that we transform O, 



comp to Ocomp \ {i} U {j}, where / G Ocomp and 



j G AO. Without loss of generality assume that / is the last column of Zj. (It cannot be an all 



zero column because, then it would not be possible for it to be in Ocomp-) Let Zi 



Z2 







Zb 


b 



Then A : 



X 

z7 



and 








b 





acting on Z achieves the transformation O, 



comp to Ocomp- 



A(Zi|Z2|/) 





T 

a 


1 


1 


^ 


ZA + bx'^ 


Zb + ba^ 





b 




x^ 


T 

a 


1 


1 





Za + bx^ 


Zb + ba^ 


b 








(z'lllz'^ll/) 



(40) 
(41) 



We claim that this same transformation can be effected by row transformations of First let 

us focus on the middle set of rows in namely the correction operators for AO. Then acting 



by M(A) 





a 





1 



gives us 



M(A)(0||0|/|Z2^) 









/ 


a 





Zl+ab^ 











1 


1 


b' 








/ 





a 


Zl+ab'^ 











1 


1 


b^' 



/|zf) 



(42) 
(43) 



Now if take the last row in Eq. ( [43] ), namely (0| |0|0| 1 = c and add xc to the top set of rows in 
Eq. ( [39l ) we obtain 



( I f ^ 



/ 







zl 



-xb' 



/ z. 



/T 



(44) 



showing the equivalence of Z and Z'. The equivalence of R and R' under A can be shown in exactly 
the same fashion as for Z\ and Z\. 
This concludes the proof that the extremal classical relations completely determine the stabilizer state. 

□ 
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5 Flipping measurement planes and temporal invariance 

In this section, we introduce a second group of symmetry transformations on MBQCs which is related 
to local complementation lITTll - |fT9l . These transformations leave the temporal relation of the measure- 
ments in a any given MBQC unchanged. An addition, if the temporal relation is a partial order, the 
distribution of the computational output is left unchanged. 



5.1 Flipping measurement planes 



According to Eq. ([3]), for any qubit a in a given resource state, the local observable measured to drive the 
computation is 



Oa[qa] 



(a) , 



-11 



a) 



Therein, is a linear function of the measurement outcomes {sb, b G Q.}, c.f. Eq. ( 15a i. 



Let's see what happens if we use a different rule for the adjustment of measurement bases, namely 



O'a [qa] = (- 1 Y'" cos (Pa O^f^ + sin % 



(45) 



That is, if = 1, to obtain C?',[l] we are flipping the observable C?a[0] about the av,^-axis rather than the 
a,|)-axis. Comparing Eqs. ([3]) and (45 1, we find that 



0'Ma\^{-iy"Oa[qal 



(46) 



independent of the measurement angle (pa- The measurements of 0'^[qa\ and are always in the 

same basis for the same qa, and the measured eigenvalues differ by a factor of (—1)^". 

We call the transformation [/] : O,- [qt] — > O'j [qi] flipping of the measurement plane at qubit i. On 
the elementary degrees of freedom, namely the resource state the Pauli observables Os'\ o^^ (action 
on ai'^ is implied), and the measurement angle (p,-, the flipping Ts\i] acts as 



(Pi 



r(0 



^ (-1)*! 



(47) 



The action of Xs\i] on the Pauli operators cr^-'' and the measurement angles (pj, for j ^ i is trivial. All our 
considerations are independent of the values of the measurement angles. In particular, the second part of 



the transformation Eq. ( 47 1 does not affect temporal order. 

We now discuss the effect of the flipping Xs[a], a G 11 on a given MBQC. Denote the measurement 



outcome of a measurement of 0'^{qa) by s'^. By Eq. (46 1, the following two measurement procedures 



are always equivalent. (I) Measuring Oa[qa] and outputting Sa, and (II) Measuring O'^^lqa] and outputting 
s[i + qa mod 2. We may call the device that performs Procedure I a (p-hox., and the device that performs 
Procedure II a s(p'-hox. Then, 





(|)-box 



s(|)'-box 
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The prime in the s(j)' -box accounts for the fact that not the measurement outcome s'^ itself is outputted, 
but rather the locally post-processed value s'^ + qa mod 2. Now, instead of outputting s'^ + qa, the device 
at a may only output s'^ (that is an s<p -box), and the classical post-processing relations for the adaption 
of measurement bases are modified accordingly, i.e. Sa — > s'^ + qa mod 2. Can the resulting relations 
again be written in a form q = T's' + H'g? 

We now attempt transforming a 0-box into an s^-box at qubit a. The vectors of measurement out- 
comes s and s' are related via s = s' -{-eae^q. Inserting this relation into Eq. ( 15a i, we obtain 



(/ + re«ej)q = rs'+//g mod 2. 



(48) 



Case I: Taa = 0. Physically, this means that the measurement basis at the flipped qubit a does not 
depend on the measurement outcome at a, before the transformation. Multiplying Eq. (48 1 with from 
the left yields qa = ejTs' + ej//g mod 2. Inserting back into Eq. (48 1, we obtain 



Likewise, 



T' = T + TeaelT mod 2. 



H' = H® TCaelH, Z' =Z® ZtaclT, R' = Re ZeatlH. 



(49) 



(50) 



Eqs. (49 1 and (50 1 completely describe the effect of the flipping Ts[a] of the measurement plane at qubit 



a on the classical processing relations Eq. ( 15a i, ( 15b). 



Remark: If the matrices H,R, T,Z are given their normal form Eq. (27 1 wrt the pair /gauge, Ocomp then 
the flipping of the measurement plane at any vertex a with Taa = leaves this normal form intact. We 



can therefore state a transformation rule equivalent to Eqs. ( 49 1, ( 50 1 for the matrices H, R, T, Z. This rule 
is, in fact, simpler. We assemble the composite matrix 



(51) 



which is a square matrix of size ( | ^1 1 + |/gauge | + 1 Ocomp I ) x ( I ^ I + 1 4auge I + 1 Ocomp | ) ■ The effect of flipping 
of the measurement plane at a then is 













I 











H 


T 








\r 


Z 


/ 


/ 



T'ext = Text + Text 6^6 J Text mod 2 . 



(52) 



Thus, the rule is just the same as Eq. (49 1 for the original influence matrix T. 

Furthermore, Text is the influence matrix for MBQC on a bigger resource state constructed from 
For the support Q.' of we require two additional sets of qubits, /' and O', with |/'| = \0'\ = 



J gauge I 



|Ocomp|. such that D,' = D,[J I' U O' . is obtained from by the following construction: 




(53) 
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Therein, the gates O — O = As are as-controlled as-gates, i.e., Ayai^^Aj = ai'\ A^a^^^Aj = a^^^^ 0aP'\ 
etc, and the extra qubits in /' and O' are initially prepared in the eigenstate of aw, with eigenvalue 1. With 



of Eq. (521 is /'I (a 



the definition Eq. ( 53 1 of , the labelling of the blocks of rows and columns for the matrix on the r.h.s. 

\0' for the rows. 



'comp 



Y\Ocomp\0' for the columns and 



1(4 



gauge I gauge 



We thus find that all information in the classical MBQC processing relations is temporal information 
for the computation on a slightly extended resource state. 

Case 2: Taa = 1. In this case, the measurement basis at a does depend on the measurement outcome 
at a. This is an example for a closed time-like curve (only involving the measurement device at a), and 
an obstacle to deterministic runnability. Now, the matrix / ( 
invertible, eJ(/( 



I Tta^a '^he left side in Eq. (48 1 is not 



0. Hence, the relation (48 1 can not be solved for q in this case. There is no 



relation Eq. ( |15a[ ) with the same sets /gauge, Ocomp before and after flipping. 

We now discuss the consequences of flipping measurement planes for the above two cases. 



5.2 Flipping measurement planes and local complementation 

We now return to the above Case 1, namely when flipping of a measurement plane yields a computation 



with a new relation Eq. ( 15a i. Note that the computation before and after the flip generate the same 
output distribution. Flipping a ^-hox into an s^'-box is an equivalence transformation, only based on the 



operator identity Eq. (46 1. Changing an s^'-box into an s0-box is again an equivalence transformation, 
provided it can be carried out. 

The influence matrices T and T' before and after the flipping, respectively, are in general not the 



same, c.f. Eq. (49 1. However, T and T' still generate the same temporal order, as we now show. 

Lemma 4. Be T an influence matrix with Ta = 0. Then, T and T' = T (B TetefT generate the same 
temporal relation under transitivity. 

Proof of Lemma^ Let's introduce a shorthand a — )• c for c G fc{a) (meaning that the measurement 
outcome at a influences the measurement basis at c). Now, we have to show that e f <^=^ ^ £^t' /, 



for any T' generated from T by the transformation Eq. (49 1. 

(I) "=^>": Assume that e f. Then there exists a sequence of measurement events e — )• mi — )• 
m2—;'..—)-a—>-c— 7". .—)•/. Can we break the arrow a — )■ c, say? To investigate this, let us rewrite the 



transformation rule Eq. (49 1 for the flip Ts[/] as 



fc(a) 
fc(a) 



fc(a) I 
fc(a), 



)fcm 



if / G fc{a), 
if i^fc{a), 



(54) 



Case 1: a — )• / before the transformation Ts[/]. Then, fc'(a) = fc(a) ©fc(/). Since Tu = by assumption, 
a — >• / after the transformation Ts[/]. Sub-case la: / — c before the transformation Ts[/]. Since Tu = 
(/ fc{i)), i — c after the transformation Ts\i]. Thus a ^ i ^ c after the transformation, and hence 
a ^T' c- Sub-case lb: / -/^ c before Ts\i]. Then, a ^ c remains after the transformation. Case 2: a -/^ i 
before the transformation Ts[/]. Then a ^ c after Ts[/]. Thus, in all cases a ^j' c, and therefore e ^ji f. 
(II) " " 



": FromEq. (54i, t4/]^ = /. □ 



Apply a series of transformations Eq. ( 49 1 on an initial influence matrix T with vanishing diagonal 
part may produce an influence matrix with a non-vanishing diagonal part. Thus, the application of the 



transformation Eq. ( 49 1 is restricted. To circumvent this problem, we introduce a modified transformation 



T' = T + TaefT + ^{TeiejT) mod 2. 



(55) 
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Clearly, this transformation takes influence matrices with vanishing diagonal part to influence matrices 
with vanishing diagonal part, and thereby avoids the problem of restricted applicability of transforma- 



tion Eq. (49 1. Note that the transformation Eq. p5h has the form of local complementation, albeit the 



influence matrix T that it acts on will in general not be symmetric. 



But what is the physical significance of transformation Eq. (55 1? The only additional effect of the 



transformation f[/] over Ts[i] is the cancelling of the diagonal part of the influence matrix after the trans- 



formation, c.f. the last term in Eq. (55 1. This can be achieved by a local unitary that exchanges a^,^ -f-)- CT, 



on a respective qubit. The action of on the elementary degrees of freedom therefore is 



■'s(l) ' 

m-- -^'r. ' ' ... (56) 



(Pi (-!)*!-<?'/, 

^ yjefc{i)nbc{i) 
\^) — > \^) 



We find that the local measured operators for all qubits j G bc{i) nfc{i) change in a way that cannot be 
accommodated by a change of the respective measurement angle. For those qubits, the new measured 



observables lie in a different equatorial plane of the Bloch sphere. Therefore, the transformation Eq. (56 1, 
unlike the transformation Eq. ( |49| ), does not necessarily map a given computation onto itself. What it 
does, however, is mapping a given computation to a computation with the same temporal relation. 

Lemma 5. Be T an influence matrix with Tu = 0. Then, T and T' = T ® TetefT® 2>{TeiejT) generate 
the same temporal relation under transitivity. 

Proof of Lemma [5] Assuming the initial influence matrix has vanishing diagonal part, we split the 
transformation T — > T ® T titj T ® {T titj T) into two steps, namely T — > T' = T ®TtitjT and 
T' — > T" = T' ® &{T'). By Lemma |4j T and T' generate the same temporal order. Now assume that 
T^f, = 1 for some k ^ Q., k ^ i. This requires that Tj^ = T^i = 1. Then, we also have 7^^ = T^'. = 1 and 
T-l = T^'- = 1. Thus, k^T' k and k -<t" k. The closed time-like curve involving k is not changed by setting 
r^'[ = 0. All other relations trivially remain unaffected by the transformation T — T' . □ 



6 MBQC - a toy model for quantum space time? 

In attempts to unify the theory of general relativity with quantum mechanics, often the viewpoint is 
taken that spacetime is not an independent construct, but rather a consequence of the laws of quantum 
mechanics. Once this assertion is spelled out, the natural next step is to identify the key quantum property 
which yields a mechanism for generating temporal order, and to illustrate this mechanism in a toy model. 

We do not solve any puzzle of quantum gravity here, but argue that measurement-based quantum 
computation possesses certain properties that one expects to find in a toy model generating spacetime 
from none. Namely, (i) Despite its origin in non-relativistic quantum mechanics, time in MBQC is a 
binary relation among spacetime events (here quantum-mechanical measurements). In computations 
without closed time-like curves — which are the ones of practical interest — this relation is a partial or- 
dering. There is no external time parameter, (ii) MBQC invokes a quantum mechanical principle that 
strongly constrains the possible temporal orders, namely that the logical processing is not affected by the 
randomness inherent in quantum mechanical measurement. This principle can be formulated in terms of 
an underlying symmetry group. 

Below we discuss aspects of MBQC that we find are of interest for toy models generating spacetime. 
These are: an analogue of Malement's theorem, forward cones arising as solutions of a wave equation. 
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and event horizons. We emphasize that these analogies arise at a formal level — MBQC is about bits, not 
matter fields. 



6.1 An MBQC counterpart of Malament's theorem? 

Malament's theorem [22J states that in any spacetime manifold the light cones determine the metric up 
to a conformal factor. Here we argue that Theorem [4] is an MBQC counterpart of that. To begin, we 
must identify an MBQC counterpart of spacetime. We say that the MBQC equivalent of a spacetime 
is the entire measurement-based quantum computation, with its resource state and temporally ordered 
measurement events. The reason for this identification is that in Malament's setting of General Relativity, 
spacetime is all there is to reason about while in our setting it is the process of MBQC. We further identify 
the spacetime points with the location of measurement events. 

Recall that an MBQC is fully specified by (i) the stabilizer generator matrix of the resource state. 



(ii) the set of measurement angles, and (iii) the linear processing relations Eq. (15a i, ( 15b I for the adap- 
tion of measurement bases and generating the output. Theorem |4] states that (iii) determines (ii), i.e. all 
of the MBQC except the measurement angles. 

We can thus read Theorem [4] as an MBQC counterpart of Malament's theorem if we make two 
further identifications: (a) the light cone structure in GR corresponds to the linear processing relations in 
MBQC, and (b) the conformal factor in the spacetime metric at every spacetime point corresponds to the 
measurement angle at every measurement location. 

Regarding (a), sure, the linear processing relations contain the forward cones in MBQC thorough 
the matrix T, but don't they also contain the additional matrices H, R and Z that have nothing to do 
with MBQC temporal order? In this regard, note that by doubling the qubits on the input and output 
boundaries /gauge and Ocomp in the support of the resource state, the influence matrix Text of the extended 



resource state comprises all information about H, R, T and Z (c.f. Section 5.1 1. In this sense, all 
information in the processing relations is temporal. 

Regarding (b), we presently do not know of a physical reason for identifying a scale factor for the 
metric at any spacetime point with a measurement angle at every MBQC spacetime point. However, we 
note that both Malamat's theorem and its MBQC analogy Theorem |4j in their respective settings, leave 
one real- valued parameter per spacetime point unfixed. 



6.2 Forward cones satisfy a wave equation 

We consider an MBQC on a resource stabilizer state in graph normal form We show that if the 
measurement plane is (X/F) for every qubit then the forward cones satisfy a wave equation. Here, 
{X /Y) is an unordered pair such that we are left with both the possibilities of = X and = Y. The 
choices may differ for different qubits. 

Following ll37l . we define the discretization of the Laplacian A = ^f^j ^ as A : = DD^ , where D is 
the incidence matrix of the graph G. In our case, all addition is mod 2. An equivalent formulation then 
is 

A:=r + ^ mod 2, (57) 

with r the adjacency matrix of G and & a diagonal matrix such that [^]vv = deg(v), for all v G V{G). 
The discretized Laplacian can act on the characteristic vector of fc{a), yielding Af^. For further use 
we define the function /« : V{G) — > Z2, by fa{v) = 1 (0) if v G fc{a) (v fc{a)). Then, the action of A 
on fa is defined through A/a(v) := [Afo],,. 
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The definition Eq. ( [57] ) of the Laplacian is in accordance with intuition, as the following example for 
two-dimensional lattice graphs (embedded on a torus) shows. There, 

{Afa){x,y) = [ifa{x+l,y)-fa{x,y))-{fa{x,y)-fa{x-l,y))] + 

+[ifa{x,y + \) -fa{x,y)) - {fa{x,y)-fa{x,y-m mod 2 
= fa{x+l,y)+fa{x-\,y)+faix,y+l)+fa{x,y-\) mod 2 
= [(r + i^)f,](,,,) mod 2. 

In the last line, ^ = because all vertices have degree 4. Back to general graphs G, the offset & in 



Eq. (58 1 is necessary such that Af = whenever / = const. We now have the following 



Lemma 6. Consider MB QC on a graph state \G) where every qubit is measured in the {X ,Y)-plane. 
Then, the forward cones satisfy a wave equation with position-dependent mass m. 



miv)faiv)+Afaiv) = 5a,„yvG0: 



comp ' 



(58) 



Therein, addition is mod 2, and m(y) = deg(y) +b(y) mod 2 with b{v) = 0(1) if the measurement plane 
at vis [X,Y] (is [Y,X]). 

Proof. For a given qubit b, irrespective of whether the measurement plane at b is [X,Y] or [Y,X], with 



Eq. ( 19 1 for the correction operations K{a) we find that b G fc{a) iff K{a)\t G Xh,Yb. In words, K{a) 
has Pauli operators ~ X exactly in those places that are in the forward cone of qubit a. Then, using the 
Definition [T] of graph states, we find 



K{a)=X{fa)Z(rfa), 



(59) 



where X(g) := (g) 

Special case: The measurement plane is [X,Y] for all qubits ^(^^"^ =Xa, Va S V{G)^. Now, for all 
by^awe require K{a)\b € {Ib,Xh}, and K{a)\a G {Za,Ya}. Using Eq. (59 1, we thus find the constraint 



(60) 



Therein, is a |C?^(,jj^p| -component vector which has an entry 1 in the ath component, and zeros every- 
where else. With Eq. (57 1, this reproduces Eq. (58 1 for m(y) = deg(v), in accordance with Lemma[6]for 
b{v) = 0. 

General case. For b ^ a we require that K{a)\i, G {Ih,Xh} if the measurement basis at b is [X,Y] 
and K{a)\i, G {Ih,Yb} if the basis at b is [Y,X]. Furthermore, we require that K{a)\a G {Za,Ya} if the 
measurement basis at a is [X,Y] and K{a)\a G {Z„,Xq} if the measurement basis at a is [Y,X]. Using 



Eq. ( 59 1 again, we find Eq. ( 58 1, with m(v) = deg(v) +b{v) mod 2 as required. □ 



6.3 CTCs of length 1 and event horizons 

In this section we consider an MBQC that has a closed time-like curce of length 1 at a given qubit /. That 
is, the measurement basis at / depends on the measurement outcome at /. We argue that if we flip the 
measurement plane at qubit / we obtain an MBQC from which the closed time-like curve is removed, but 
instead qubit / vanishes behind the MBQC counterpart of an event horizon. 

Be / a qubit such that Tu = 1 before flipping the measurement plane at qubit /. Such a qubit / cannot 
be in /gauge, since /gauge ^ / by definition. If Tu = 1 then / G bc{i). The backward cones of all qubits in / 
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are empty by definition of /, liowever. Likewise, / ^ Ocomp- If Ta = 1 then / e fc{i). However, fc{a) = 
for all a G Ocomp- Thus, there is only one case to consider, namely / G {LaugeY H (Ocomp)'^- 

In this case, there exists a correction operator K{i) for qubit / before the flipping, K{i) = a^^ (g) 
K{i)\£i\i- After flipping at /, this operator turns into 

TMm) = of ®m\^\, =: r (/). (6i) 



That is, the operator {K{i)) resulting from flipping at / is a gauge type operator, c.f. Eq. (22 1. Thus, 
the flipping transformation (when Tu = 1) enlarges /gauge by one qubit, 

Tv[''] • /gauge ^ /gauge U if Ta — 1. 

Furthermore, after the flipping at / there no longer is a correction operation for qubit /, hence 

: Ocomp — > Ocomp U {/}, if Tu = 1. 

This has two consequences. First, the forward cone of / becomes empty. In particular Tn = after the 
flipping. Thus, the closed time-like curve consisting of qubit / has been removed. Second, an additional 
bit of optimal classical output is being created by the flipping at /. What does that output bit signify? 
Recall that before the flipping at /, the rule for adjusting the measurement basis at / is 

! 

1i = Si+ mod 2, 

for some set J (^Q.. Here, we have dropped a constant offset h^g on the r.h.s. The symbol "!" above 
the equality means that equality is a requirement for the correctness of the computation, but it cannot be 
deterministically implemented. As follows from Eq. (|46]), the measurement outcomes before and after 
the flip. Si and 5- are related via = ^i®qi. For all the other qubits, s'j = sj. Substituting this into the 
above relation, we obtain 

I 

s'i+ 4 mod 2 = 0, V^7^gZ2. (62) 

Thus, the additional output bit o,- = s\ + Y.jeJ\i^'j mod 2 is a flag. If = then the computation suc- 
ceeded, and if o, = 1 then it did not. 

Now suppose that the problem solved by the given MBQC is in NP. Then, this flag bit is not neces- 
sary. The remaining output may be efficiently checked for correctness anyway. Thus, one may safely 
discard the extra bit of output. Not post-selecting on o,- = can, if anything, only increase the success 
probability of the computation. We thus arrive at 

Lemma 7. Be an MBQC with a classical output o and influence matrix T such that Tu = 1, i.e., .M\ 
has a closed time-like curve involving a single qubit i G D.. Be .M2 the MBQC with the same classical 
output o, obtained from ^\ by flipping the measurement plane at i. Then, the closed time-like curve 
of i in ^\ is removed in Furthermore, if ^\ solves a problem in the complexity class NP with 
probability p, then ^2 solves the same problem with probability > p. 

Remark: Lemma|7]does not guard against the inefficiencies of post-selection, in particular if multiple 
CTCs of length 1 are being removed. While the success probability after removing the CTCs is guar- 
anteed not to be smaller than for the original computation with the CTCs (which can only be executed 
using post-selection), neither it is provably significantly larger. 
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Figure 3: Breaking a closed time-like curve of length 1. The looped qubit becomes an element of 
4auge ^ ^comp after flipping the measurement plane. As such both its forward cone (Ocomp ^ O), and 
backward cone (/gauge ^ /) must be empty. 



Event horizons. Let us consider the flow of information between qubit / whose measurement plane 
has been flipped and the other qubits. Before the flip (MBQC of Lemma |7j), / E (Ocomp)'^ H (/gauge )'^- 
After the the flip (MBQC of Lemma|7]), / G Ocomp H /gauge- In ^2, since / € /gauge> no information for 
the adaption of measurement basis is flowing into site / from the other sites. Likewise, since / G C?comp> 
no information for the adaption of measurement bases is flowing out of site /. Finally, because of the 
normal form Eq. ([27]), the measurement outcome 5, appears in only one readout bit. This readout bit 



is Oi as given in l.h.s. of Eq. (62i, which is precisely the bit of classical output that can be discarded 
if the problem solved by the quantum computation is in NP. If o, is discarded, then no information is 
flowing out of the site / at all. Thus, in summary, from the viewpoint of classical processing, qubit / in 
^2 becomes entirely disconnected from the computation. It vanishes behind the MBQC counterpart of 
an event horizon. 



7 Conclusions and outlook 

In this paper we have studied the constraints on temporal order in measurement-based quantum compu- 
tation which arise from the principle that the randomness inherent in quantum measurement should not 
affect the logical processing. We have established a classification of temporal relations consistent with 
a given resource stabilizer state and set of measurement planes. Conversely, we have shown that the 
linear processing relations in measurement based quantum computation, subject to the above principle, 
specify the resource state and set of measurement planes up to equivalence. We identified gauge degrees 
of freedom which need to be included in order to establish the above results. Furthermore, we found a 
transformation that leaves the temporal order in every MBQC invariant and is related to local comple- 
mentation. Finally, we pointed out formal MBQC analogues of a result and a piece of phenomenology 
in the theory of General Relativity, namely of Malament's theorem and event horizons. 
At this point, we are led to ask the following questions: 

1. We introduced a group of gauge transformations Eq. ([TSj), and a group symmetry transformations 



Eq. (55 1, generated by flipping measurement planes. Both transformations preserve MBQC tem- 



poral orders. Can the two groups be unified? 

2. Some of the temporal relations admitted by the matroid contain closed time-Uke curves. 

Given a stabilizer state and set of measurement planes Z, can we find an algebraic (or other) 
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structure which comprises only the partial orders? Can the partial order of measurements with the 
smallest set Ocomp be efficiently computed? 

3. The generators of the resource state stabilizer commute by definition. As pointed out in the second 
remark below Theorem[4j this commutativity condition constrains the MBQC linear processing re- 
lations, including temporal order. What is the physical meaning of this constraint to the processing 
relations? 

4. In MBQC, the link between the randomness in quantum mechanical measurement and temporal 
order is the principle that the randomness of measurement outcomes should not affect the logical 
processing. In a context more general than quantum computation, what could this principle be 
replaced by? 
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A Invariance under the gauge transformations 

A.l Gauge transformations and temporal order 

In this section we provide a different angle at Theorem [T] namely we show that the gauge transformations 
Eq. ( 1 8 1 impose severe constraints on the possible temporal orders for a given stabilizer state and set of 
measurement planes. 
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Symmetry constraints on temporal order in measurement-based quantum computation 



The gauge transformations Eq. (fTSll act on q and s, and can therefore have a non-trivial effect on the 



classical processing relation Eq. ( |15a ), q = Ts+Hg. We now study this effect. Since the transformations 
Eq. (18 1 are caused by the insertion of stabilizer operators into the state overlap in Eq. (|5]), they do not 
change the physical situation. Therefore, the temporal relations before and after any such transformation 
must be equally valid, although not necessarily identical. By insertion of the stabilizer operator into 
the overlap (Oiod*!'), the stabilizer of and the the sets /gauge > Ocomp do not change. Therefore, by 
Theorem[T| the matrices T and H do not change. Thus, besides q and s, all that can change in the relation 
Eq. ( 15a| ) under a transformation Eq. (18 1 is g. The following two viewpoints are always equivalent: 
(A) The relation q = /g(s), under the action Eq. ( 18 1 of a Gk on (s,q) is changed into an equivalent 
such relation q = /g'(s), with Gk '■ g — > g'- (B) The relation q = /(s,g) remains invariant under all 
transformations Gk, acting on the triple (s,q,g). We choose the latter viewpoint. 

We now infer the action of the transformations Gk on g. Without loss of generality we assume that 
the relations Eq. (|15a|) are given with H in its normal form Eq. (pT]), 



H 



g mod 2. 



Furthermore, we assume that the pair /gauge, Ocomp is extremal. Since /gauge ^ / by definition, for each 
i G /gauge> Qi only depends on g but not on the measurement outcomes s, qt = gi. Now, the correction 
operators K{a), a £ (Ocomp)'^ derived from the normal form Eq. (26) of K{a)\j^,^^^^ has no a^-part. 



Therefore, the corresponding transformations Gk do not flip qt, for all / € /gauge- In order to preserve the 
relation qi = gi, they thus leave g unchanged. Now consider the other stabilizer generators, K{i), i G /gauge, 
obeying the conditions Eq. (22 1. By construction, flips qt but no other qj, for / / 7 e /gauge- Hence, 
to preserve the relations qi = gi, it must also flip gi, but no other gj, i j € /gauge- Thus, for a stabilizer 



element K = ^aeci{<^^"^y"{<^, 



Gk: g 



■ g© w / 



(63) 



Therein, we have assumed that the basis choice for g is such that the matrix H appearing in Eq. (15a I is 
of normal form Eq. (p7|). 



We have now fully specified the action of Gk on the triple (q,s,g), c.f. Eq. (18 1, (63 1. MBQCs 
satisfy the invariance condition 



q = rs + //g mod2^^GK{q) = TGKis)+HGKig) mod 2, V/^ G 



(64) 



It is evident that the requirement (64i of invariance of the processing relations ( 15a 1 under the gauge 



transformations poses constraints on the possible matrices T and H. In fact, as we show below, given 
Ocomp the matrices T and H are uniquely specified uniquely by the above invariance condition. 

To check the invariance condition Eq. ( [64| ) in a specific case, we return to our 3-qubit cluster state 
example of Section |2.2[ We consider the effect of the transformations induced by generators Ki 



(T^^ci^^ K2 = ai^^'a^aP^ and ^3 = ai^'a^^^ on the processing relations Eq. ( |l6| ). As noted earlier, 
{1}. Then, with Eqs. (fTsj) and (|63|l. 



J gauge 



Gki 

Gki 
Gk3 



q^q©(l,0,0)^ s 
q^q©(0,l,0)^, s 
q^q©(0,0,l)^, s 



i(0,l,0)^ 
'(1,0,1)^ 
'(0,1,0)^ 



^1' 
^1- 



(65) 
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It is easily checked that the relation Eq. ( 16 1 is invariant under the transformations G/f, , G^, and Gk^ of 
Eq. ( 65 1. However, if the transformations are restricted to q, s, the relation Eq. ( 16 1 is no longer invariant 
under the transformation induced by Ki . 



We now return to the general case and show that, given the set Ocomp and the action Eq. ( 18 1, (63 1 of 
the gauge transformations on the triple (q,s,g), the invariance condition Eq. (64i uniquely specifies the 
classical processing relations Eq. ( 15al for the adaption of measurement bases. 

Recall that we write the stabilizer generator matrix for in the 00 /as-basis as = (<I>| \S). 

Then, for the stabilizer generator £ corresponding to the a-th row of with Eq. ( 18 1 

the action of the gauge transformation Gk„ on s, q is 



Gk„: s — >seroWa(5), q — ^ q row^ (O) . 



With Eq. (63 1, the action of Gk„ on g is 



(66) 



(67) 



Here, roWa(<I>)|/^^^„^ denotes row„(<I>) restricted to /„, 
q = Ts + Hg under Gk„ becomes 



Then, the condition Eq. ( 64 1 for invariance of 



roWa(<I>) = TmWa{S) +HicoWa{^)\i^.^^^^ mod 2. 
This condition must hold for all stabilizer generators Ka simultaneously, hence 

= TS^ +H^^\i „ mod 2. 



(68) 



By definition, the qubits in /gauge have empty backward cones, and the qubits in Ocomp have empty 
forward cones, hence T is of the form 








T 






where the column split is (Ocomp)'^ |Ocomp and the row split is /gauge I (/gauge )'^» c.f. Eq. (29 1. By right- 
multiplication of relation Eq. (68 1 with a suitable matrix, we transform <I*^|/„;,„„exn irito a matrix of form 
(/|0) where the column split is between /gauge and /gauge- By definition of /gauge> such a transformation is 
always possible. Under the same transformation. 













' '^^l(Ocomp)^xn > {Sl\S2). 



(69) 



Inserting the above into Eq. (68 1, we find that H must be of normal form Eq. {27), H 



, and 



<I>i = TSi+H mod 2, 
<I>2 = TS2 mod 2. 



(70) 



Now, S2 must be an invertible matrix. This is the condition that, by definition of Ocomp, every measure- 
ment outcome in (Ocomp)'^ is correctable. Then, 



T = <I>252"^ mod2, H = <I>i+<I>2>S'2"'5i mod 2. 



-1 , 



(71) 



Hence the relation q = Ts + Hg is uniquely specified. 
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A.2 Gauge transformations and computational output 



In addition to Eq. 64 we also require invariance of the classical output under the transformations Eq. ( 1 8 1, 
o = Zs + Rg mod2 = ZGK{s)+RGK{g) mod 2, V/T G (72) 



(|63ll, 



Like Eq. ( 64 1, Eq. ( 72 1 is a determinism constraint. If for a single output bit o the relation o = z^s + r^g 



is not invariant under all gauge transformations Eq. ( [18] ), ( [53] ), then the value of o is guaranteed to be 
random, and thus useless as readout bit of a computation. Specifically, 

Lemma 8. Assume an MBQC where the relation q = Ts + Hg is invariant under the gauge transfor- 
mations Eq. ( 18 1, {63 ), but an output bit o exists whose defining relation o = z^s + r^g is not invariant 
under the action of GKfor some K € S^{\^)). Then, the value of o is completely random, independent 
of the choice of measurement angles. 

Proof of Lemma^ For simplicity, consider first the special case where K e acts trivially 

on g, G/s:(g) = g, Vg. We may then write o = Y^ieJ^i + for ™ offset c = r^g. We call the string 
s|7 of measurement outcomes on J even (odd) if it has even (odd) weight. We denote the local post- 
measurement state on qubit a by | (p^ , 5^ , (s , g) ) , where (pa is the measurement angle, Sa the measurement 
outcome and qa specifies the chosen measurement basis. 

Under the transformation Gk, s — > s© As/f, where, by assumption, Asj^jy is odd. Now, the proba- 
bility of outputting o = c is 



p{o = c) 




a{%,Sa,qais,g 




a{%,Sa® Ask, a , qa © AqK,a 



a{(pa,Sa®AsK,a,qa{s®ASK,g 




(73) 



a{%,Sa,qa{s,g)\ l"^) 



Thus, p{o = c) = p{o = c) = 1 /2. Note that in transitioning from the third to the fourth line of Eq. (73 1 
we have used the invariance property Eq. (64i, i.e., the assumption that the adaption of measurement 
bases is deterministic. 

In the general case, Gk '■ s — > s © As a-, g 

1 



g, and thus p{o = c) 



2 1 'gauge 



I I 

g s|j=even 



the above we then find p{o = 0) = p{o = 1) = 1/2. □ 




g © Ag^. We note that we can choose any gauge fixing 

\ 2 



a{(Pa,Sa,qa{s,g)\ 



aeO. 



By an argument analogous to 



